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Centrality and pseudorapidity dependence of elliptic flow for charged hadrons in 
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This paper describes the measurement of elliptic flow for charged particles in Au+Au collisions at 
ys^7 = 200 GeV using the PHOBOS detector at the Relativistic Heavy Ion Collider (RHIC). The 
measured azimuthal anisotropy is presented over a wide range of pseudorapidity for three broad 
collision centrality classes for the first time at this energy. Two distinct methods of extracting the 
flow signal were used in order to reduce systematic uncertainties. The elliptic flow falls sharply with 
increasing \rj\ at 200 GeV for all the centralities studied, as observed for minimum-bias collisions at 
^/sTT = 130 GeV. 
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It is widely accepted that a very dense and pos- 
sibly new state of matter is being created in central 
Au+Au collisions Q at the Relativistic Heavy-Ion Col- 
lider (RHIC) at Brookhaven National Laboratory. The 
azimuthal anisotropy in the distribution of produced 
particles ("flow") is a consequence of the initial spatial 
asymmetry of the collision zone and subsequent rescat- 
tering processes which convert this to a final momentum 
anisotropy. Measurements of flow are therefore sensitive 
to the system early in the collision and to its dynamical 
evolution. 

There is an extensive data set on flow results from 
RHIC USS HIE II II HE EH!!, but as of yet, the 
least understood result is the pseudorapidity and energy 
dependence of the ellip tic flow, V2 (??) , measured over an 
extended grange [11L Recently some theoretical 

progress has been made by introducing a longitudinal 
dependence in the source shape and/or by assu ming 
incomplete thermalization away from r\ = [l3L Ilif . 
In this paper, we extend our earlier measurements by 
examining the dependence of the v^ij]) shape on the 
collision centrality. The elliptic flow of charged hadrons 
has been studied using data from the PHOBOS detector 
during the 2001 Au+Au run of RHIC. In addition to the 
"hit-based" method previously used in 0] , a new "track- 
based" method was developed and employed to improve 



the accuracy of the measurement and provide a valuable 
consistency check of the hit-based analysis. 

The PHOBOS detector consists of silicon pad detec- 
tors arranged in single and multiple- layer configurations 
surrounding the interaction region, as described in [i"o| . 
The two multiple-layer magnetic spectrometer tracking 
arms are configured with a field free region near the 
interaction vertex followed by tracking inside the magnet. 
This leads to two classes of found tracks with different 
acceptances. "Straight-line" tracks cover < rj < 1.8 
with an azimuthal acceptance of A(j> w 22° centered at 
= 0° and 180°. "Curved" tracks cover < r) < 1.5 
with a variable azimuthal acceptance of A<j> ss 20°. 
Details of the tracking procedure are given elsewhere . 
The single layer configuration includes the octagonal 
multiplicity detector (OCT) with |?7| < 3.2 and six an- 
nular silicon ring multiplicity detectors (RINGS), with 
3.0 < \i]\ < 5.4. The rings and most of the octagon have 
full azimuthal coverage except near the middle of the 
detector (mid-rapidity for nominal vertices) where the 
azimuthal coverage drops by a factor of two. 

Two sets of scintillating paddle counters were used for 
triggering and centrality determination 0, 0, 0] . In 
addition an online vertex trigger was employed, using 
two sets of Cerenkov detectors. The hit-based method 
required events whose collision vertex (v z ) was centered 
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FIG. 1: Elliptic flow as a function of pseudorapidity (t>2(f?)) 
for charged hadrons in minimum-bias collisions at ^/s NN — 
130 GeV (open triangles) ^3 an d 200 GeV (closed triangles). 
One sigma statistical errors are shown as the error bars. 
Systematic errors (90% C.L.) are shown as gray boxes only 
for the 200 GeV data. 



at —34 cm away from the nominal vertex position, 
along the beam axis The vertex trigger enabled 

a special sample (~ 1 million triggers) of such events 
to be taken. The track-based method required events 
with vertices within about 10 cm of the nominal vertex 
position, which allowed a large fraction of the 2001 
Au+Au data set at 200 GeV (~ 25 million triggers) to 
be used. The minimum-bias sample for the hit-based 
method consists of all triggered events that have a valid 
reconstructed vertex. This engenders biases similar to 
those discussed in [ijj and leads to the average number 
of participants (N part } given in Table [I] For the track- 
based method, only the fraction of the cross-section 
unbiased by trigger and vertex inefficiencies is used to 
form the minimum-biased sample. The average number 
of participants for this method is also given in Table Q] 
For the centrality dependent Vzirf) analysis, the data 
samples were subdivided into the three centrality classes 
given in Table [I] The top 3% of the cross-section, 
where the flow signal is smallest, was omitted to reduce 
the resulting statistical and systematic errors on the 
most central bin. Differences in the average number 
of participants between the two methods, for the same 
fraction of the Au+Au cross-section, occur because the 
track-based method is track weighted whereas the hit- 
based method is event weighted. This results in slightly 
higher (N part ) values for the track-based method, which 
are insignificant given the systematic error in (N part ). 
For both methods the resulting centrality classes are 
unbiased. The summary of the number of events used 
is also given in Table [I] 

Monte Carlo (MC) simulations of the detector per- 
formance based on the Hijing [2(| event generator and 
GEANT 3.2lj23 simulation package were used for sys- 
tematic error studies. 

Figure ^ shows the minimum-bias result for the 



200-GeV data using the hit-based method (as described 
in 23). The data show a steady decrease in w 2 with 
increasing \rj\, similar to that seen at the lower energy of 
y/s NN = 130 GeV (also shown). No significant difference 
in shape or magnitude is seen within the systematic 
errors. The ratio of V2 at y/s NN = 200 GeV compared 
to 130 GeV, averaged over all 77, is 1.04+ 0.03(stat.) ± 
0.04(syst.). 

The track-based method correlates the azimuthal angle 
of tracks that traverse the spectrometer, 4>t r ki with the 
event plane as measured in the octagon, ^2, event by 
event. The method used is based upon the scheme 
described by Poskanzer and Voloshin (2^, where the 
strength of the flow is given by the n th Fourier coefficient 
of the particle azimuthal angle distribution 
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l + ^2w„cos[n(^ rfe -* R )]. (1) 



In this analysis only the n = 2 component is studied and 
the true reaction plane, ^> r, is approximated by the event 
plane 'J' 2- 

The use of tracking requires events with vertices 
near the nominal vertex range (—8 cm < v z < 10 cm) to 
ensure maximum track acceptance in the spectrometer. 
Only the parts of the OCT detector with complete 
azimuthal acceptance (i.e. those away from mid-rapidity) 
are used to determine the reaction plane. Two sub- 
events, symmetric in rj and of equal charged particle 
multiplicity, are used to determine the event plane 
resolution. The sub-event sizes are vertex dependent, 
resulting in a resolution correction that is both centrality 
and vertex dependent. The resulting sub-event ranges lie 
between 2.05 < |7y| < 3.2, and are widely separated, thus 
greatly reducing the effects of any short range non-flow 
correlations. The event plane is determined using 
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J2i u>i sin(2<fc) 
Y.iWi cos(2(f>i) 



(2) 



where <f>i is the i th hit's measured angle, and the sums 
run over all hits in both sub-events. The sub-events are 
combined for the event plane determination in order to 
maximize its resolution. Vertex dependent corrections, 
some determined on an event-by-event basis, are used 
as weights (wi) ^3 m order to remove acceptance and 
occupancy biases. The resulting distributions of event 
plane angles are found to be flat within 2%. 

To determine the V2 coefficient, the measured 
distribution is divided by a mixed event 
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distribution in order to remove detector related effects, 
such as non-uniformities in the azimuthal acceptance of 
the spectrometer: 
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where A(f> denotes 
resolution correction. The 



~l + 2(-^)cos(2A0), 
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^>2 and Cres is the event plane 
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distribution, with 
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Hit-based 




Track-based 




/0<J Au+Au 


< Npart > 


Number Events 


/OU Au+Au 


< Np ar t > 


Number Events 


minimum-bias 




205 


34,727 


0-50 


236 


5,050,778 


central 


3-15 


288 


11,221 


3-15 


294 


1,439,923 


mid-central 


15-25 


199 


7,550 


15-25 


202 


1,230,394 


peripheral 


25-50 


111 


10,127 


25-50 


115 


3,087,599 



TABLE I: Characteristics of the event samples used in the two flow analyses of 200-GeV Au+Au collisions. The systematic 
error in {N par t) is approximately ±4 participants. 



zero flow, is constructed using an event mixing technique, 
where the <f>track of tracks in one event are subtracted 
from the \&2 of another event. 



Normalized A<p distributions, j^- 



mcasurcd i 



dN 
dAd, 



are determined as a function of vertex position 
and for fine centrality bins (~ 5% of the cross-section per 
bin) since the event mixing technique requires similarity 
of the class of events examined. The centrality bin 
and vertex dependent event plane resolution correction 
CVes (centrality, v z ) are determined using the sub-event 
technique |22| as 

1 

C r es (centrality, v z ) 



V2c 



cos 



Or - ^ >o 



(4) 

where if)^ ancl lt p2 > are the event planes from each 
sub-event. The ay/2 factor converts the single sub- 
event resolution correction into a combined sub-event 
resolution correction. The a factor is sometimes approx- 
imated as unity, but this approximation can break down, 
particularly when the event resolution is good. Its exact 
form is given in Refs. [22j and |23j. For the resolutions 
measured in this data set, 0.95 < a < 1. 

After averaging Eq. [3] over vertex positions and cen- 
tralizes falling into each broad centrality class defined 
in Table m the v-i coefficient is extracted from the fit to 
an even-harmonic series. (Including orders higher than 
n = 2 in the fit did not effect the extracted v%-) It should 
be noted that the resulting vi is a track weighted result 
over the broad centrality classes, since limited statistics 
precluded the Vi from being determined for each fine 
centrality bin and then event weighted. For further 
details on this technique see |23j. 

Extensive MC simulations have shown that the magni- 
tude and shape of the flow signal are correctly reproduced 
by this method. No further corrections to the measured 
i>2 coefficient are necessary, such as potential corrections 
due to the density of particles or suppression corrections 
due to backgrounds, as required in the hit-based method. 

In addition to the sources of systematic errors con- 
sidered for the hit-based analysis 11], other studies 
performed for the track-based method include analysis 
of the effects related to tracking, such as varying cuts 
on the distance of closest approach of tracks to the 
collision vertex, differences between results obtained from 
the two spectrometer arms, momentum resolution, and 
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FIG. 2: Elliptic flow (uafljjl < 1)) as a function of {N par t) 
determined by the track-based method (closed circles) and 
hit-based method (closed triangles) for Au+Au collisions at 
200 GeV. The open triangles are the results from Au+Au 
collisions at 130 GeV. One sigma statistical errors are shown 
as the error bars (within the points for the track-based 
method); gray and open boxes show systematic uncertainties 
(90% C.L.) for the 200-GeV results from the hit-based 
and track-based methods, respectively. The line shows a 
calculation from hydrodynamics |24| at x/s NN = 200 GeV. 



dependence on the bending direction. Additionally, con- 
tributions due to the vertex dependency of the resolution 
corrections, different beam orbit conditions, and errors of 
the fit parameters were also accounted for. 

Figure El shows the centrality dependence of the «2 
determined using the straight line track-based method 
over a range of < rj < 1, allowing a direct comparison 
with the same result using the hit-based technique. The 
two techniques agree well over the full range of centrality. 
The curve in Figure |21 shows a hydrodynamic calculation 
for Au+Au collisions at ^/s^ = 200 GeV . As seen 
for Au+Au collisions at 130 GeV Q (open triangles), 
the 200-GeV results at mid-rapidity are consistent with 
expectations from hydrodynamic models. There is no 
significant difference between the 130- and 200-GeV data 
in either the shape or magnitude of v-i at mid-rapidity as 
a function of centrality within the errors. 

Using tracks that traverse the full field region of the 
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FIG. 3: Elliptic flow as a function of transverse momentum 
(v2{pt)) for charged hadrons with < r\ < 1.5 for the most 
central 50% of the 200 GeV Au+Au inelastic cross section. 
The one sigma statistical errors are shown as the error bars. 
The gray boxes represent the systematic errors (90% C.L.). 
The data points are located at the average pr position within 
a pr bin whose size is given by the horizontal error bars. The 
curve shows a calculation from hydrodynamics |24^|. 



spectrometer, the transverse momentum dependence of 
the flow strength V2{pt) can be measured. This is 
shown in Figure for the top — 50% centrality for 
tracks averaged over the range < i] < 1.5. The curve 
shows the prediction of a hydrodynamical model [24| . As 
previously observed 0, the v 2 rises as pr increases and 
at pt above 1.5 GeV/c tends to flatten out well below 
the hydrodynamic curve. 

In these analyses, the reaction plane is determined 
in sub-events that are at different pseudorapidities from 
those where the V2 is measured. This should significantly 
reduce the contribution of any non-flow effects to the 
measured 1*2, particularly those due to short-range corre- 
lations. Comparisons of the V2(pt) result to the reaction 
plane and cumulant methods results from reference 
averaged over a similar centrality range, show that our 
result is most consistent with the one obtained with the 
four particle cumulant method |2^, suggesting that our 
track-based methodology is indeed largely immune to 
non-flow effects over the range |?7| < 1.5. 

Figure 0] shows «2(?7) for three centrality classes as 
defined in Table [I] Excellent agreement is seen across 
all of the centrality classes over the range of overlap 
suggesting that our hit-based method is also minimally 
affected by non-flow effects around mid-rapidity. 

To examine how the shape of the distribution changes 
with centrality, the results of the hit-based method and 
track-based methods are combined. Although obtained 
in the same experiment, the measurements should effec- 
tively be considered independent of each other due to 
the very different methods and elements of the PHOBOS 
detector used; hence the results for each method are 
combined with the reasonable assumption that the errors 
are uncorrelated. First, the hit-based results that are ap- 
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FIG. 4: Elliptic flow as a function of pseudorapidity (112(77)) 
for charged hadrons from 200- GeV Au+Au collisions for the 
three different centrality classes described in the text, ranging 
from peripheral to central (25 - 50%, 15 - 25%, 3 - 15%) 
from top to bottom. The triangles are the results from 
the hit-based method, and the circles are from the track- 
based method. The open circles are the track-based results 
reflected about mid-rapidity. One sigma statistical errors are 
shown as the error bars (within the points for the track- 
based method) ; the gray and open boxes show the systematic 
uncertainties (90% C.L.) for the hit-based and track-based 
methods, respectively. 



proximately an equal rj distance away from mid-rapidity 
are combined (e.g. 77 = —4.87 with 77 = +5.06), weighted 
by their statistical uncertainties. The points at (77) = 
—3.05 and (77) = +3.67 are just reflected due to the lack of 
symmetry of these points around 77 = 0. The track-based 
results at 7/ = +0.17 and 77 = +0.46 are also combined 
to give V2 at 77 = +0.31, and similarly for 7/ = +1.05 
and 77 = +1.34, to give a at 77 = +1.20. The hit-based 
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FIG. 5: Elliptic flow as a function of pseudorapidity («2(??)) 
from 200-GeV Au+Au collisions for the three centrality bins 
(3— 15% circles, 15 — 25% triangles, 25 — 50% squares). Data 
for r) > are determined by reflecting the hit-based results 
about mid-rapidity and then combining them with the track- 
based results and are shown with the corresponding combined 
90% C.L. statistical and systematic uncertainties. The same 
data are reflected around r\ = and shown as open symbols. 
In the range where the methods overlap, the insert shows the 
ratio of the peripheral to central results, with the appropriate 
90% C.L. combined uncertainties. 



results and the track-based results with similar r\ binning 
are averaged, weighted by their combined statistical and 



systematic uncertainties. The resulting data are shown 
in Figure^ The pseudorapidity dependence of v-i for the 
3 centrality bins is similar to that observed in Fig^for 
minimum-bias data. For peripheral collisions, v-i clearly 
already has a non-zero slope over the range — 2 < rj < 2. 
The overall shape of ^2(77) is not strongly centrality 
dependent within the uncertainties, appearing to differ 
only by a scale factor. This is illustrated in the insert 
of Fig. which shows that the ratio of the peripheral 
to central data around mid-rapidity is approximately 
constant. However, it should be noted that the central 
data around mid-rapidity is also consistent with a flat 
distribution, given the uncertainties. 

In summary, we have measured the centrality depen- 
dence of V2(rj) in Au+Au collisions at y/s NN — 200 GeV. 
Excellent agreement with the track-based method further 
validates the use of the hit-based method. This method 
allowed for the study of the i>2(??) dependence over the 
large range of 77 covered by the PHOBOS single-layer 
silicon detectors. The 200-GcV results clearly show that 
vi decreases with increasing |^|, as seen for the 130- 
GeV Au+Au collisions. From comparisons of the vi{j>t) 
results with four particle cumulant results we conclude 
that our flow measurements are largely immune to non- 
flow effects, over the range |?/| < 1.5. 

The predominant features of the viiif) distribution do 
not change significantly as a function of centrality from 
(N par t) ~ 290 to {N part ) ~ 110. The flow still falls off as 
one moves away from mid-rapidity. It is hoped that this 
data can be used to more fully understand the strong 77 
dependence of the 1*2 flow component. 
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